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ABSTRACT 

Thermal and chemical evolution of primordial gas clouds irradiated with far- 
ultraviolet (FUV; hv < 13.6 eV) radiation is investigated. In clouds irradiated by 
intense FUV radiation, sufficient hydrogen molecules to be important for cooling are 
never formed. However, even without molecular hydrogen, if the clouds are massive 
enough, they start collapsing via atomic hydrogen line cooling. Such clouds continue 
to collapse almost isothermally owing to successive cooling by H~ free-bound emission 
up to the number density of 10 16 cm~ 3 . Inside the clouds, the Jeans mass eventually 
falls well below a solar mass. This indicates that hydrogen molecules are dispensable 
for low-mass primordial star formation, provided fragmentation of the clouds occurs at 
sufficiently high density. 

Subject headings: cosmology: theory — early universe — galaxy formation — molecular 
processes — stars: formation 

1. Introduction 

Since the dawn of evolutionary cosmology, great efforts toward understanding the formation 
of the first structures have been made by many authors (e.g., Saslaw & Zipoy 1967; Peebles & 
Dicke 1968; Matsuda, Sato, & Takeda 1969; Doroshkevich & Kolesnik 1976; Carlberg 1981; Palla, 
Salpeter, & Stahler 1983; Kashlinsky & Rees 1983). In particular, partly stimulated by observa- 
tional advances and by the consequent detectability of the first stars by next-generation facilities, 
this area is being vigorously investigated these days (e.g., Haiman, Thoul, & Loeb 1996; Ostriker 
& Gnedin 1996; Tegmark et al. 1997; Yamada & Nishi 1998; Bromm, Coppi, &i Larson 1999; Abel, 
Bryan, & Norman 2000; Fuller & Couchman 2000). Those studies have repeatedly emphasized the 
key role played by molecular hydrogen. 

The role of H2 in the cosmological context is twofold. First, small primordial clouds whose 
virial temperature is 1000 < T v \ v < 8000K can cool and start collapsing via H2 line cooling after 
virialization (e.g., Haiman et al. 1996). Second, the Jeans mass inside a collapsing primordial cloud 
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is eventually reduced to much less than a solar mass as a result of H2 cooling (e.g., Palla et al. 
1983). 

Regarding the former role, there is some dispute over the role of H 2 in the formation of second- 
generation stars, namely, what happens after the first stars have formed by H2 cooling. Formed 
first stars or first quasars can be sources of UV radiation. Whereas ionizing photons are trapped in 
an HII region around the source, far-ultraviolet (FUV; hu < 13.6 eV) photons 1 travel farther and 
form an outer H2 photodissociation region (PDR). Therefore, before the overlap of the HII regions, 
or in other words the reionization of the universe, the H2 PDRs must cover the whole intergalactic 
space. At that time, the universe is filled with FUV background radiation (e.g., Haiman, Rees, & 
Loeb 1997; Ciardi, Ferrara, & Abel 2000). 

Hydrogen molecules in small objects can easily be photodissociated by either internal or exter- 
nal UV sources. Haiman et al. (1997) found that in the presence of UV background radiation at the 
level needed to reionize the universe, molecular hydrogen is photodissociated in small cosmological 
objects. Furthermore, Omukai & Nishi (1999) pointed out that even a single O star formed in a 
small cosmological object is able to photodissociate the whole original cloud. In addition, those 
clouds are easily blown out by a few supernovae because of their shallow potential wells (e.g., Mac 
Low & Ferrara 1999). These effects strongly regulate the subsequent star formation within small 
pregalactic clouds. 

In contrast to this, the temperature of sufficiently massive clouds reaches above ~8000K, where 
atomic Lya emission begins to be an effective cooling agent. Regardless of the H2 photodissociation, 
such clouds can start collapsing dynamically. It is frequently assumed that, even in primordial 
clouds that cool by the atomic cooling instead of the H2 one, star formation occurs in the same 
way as in those with the H2 cooling. However, it is far from evident. 

Recall another role of hydrogen molecules in primordial star formation. The formation of pro- 
tostars in primordial clouds with H2 cooling has been investigated by many authors (e.g., Matsuda 
et al. 1969; Yoneyama 1972; Carlberg 1981; Palla et al. 1983; Omukai & Nishi 1998). It is known 
that those clouds cool by H2 rovibrational line emission (for number density n < 10 14 cm~ 3 ), H2 
collision- induced emission (for 10 14 < n < 10 16 cm~ 3 ), and the H2 dissociation (for n < 10 20 cm~ 3 ). 
Eventually, a very small protostar (~ 1O~ 3 M ) forms at the center when n ~ 10 22 cm -3 . 

On the other hand, we know little about the later evolution of clouds that start collapsing via 
atomic cooling. It is still unclear under what condition H2 eventually forms, or whether it plays 
some role as a coolant in such clouds. It is also not known how small a Jeans mass is reached 
through atomic cooling alone. If this value falls well beyond the mass of normal stars (for example 
if it corresponds to the size of super massive stars or very massive objects), we must conclude 
that normal stars does not form only by atomic cooling. The aim of this paper is to answer these 
questions. 



1 Some authors use the term "soft-UV", instead of "far-UV", for radiation whose energy is below the Lyman limit. 
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In this paper, we study the thermal and chemical evolution of primordial clouds irradiated 
with FUV radiation. Thereby, we demonstrate that the Jeans mass inside the clouds is reduced 
to much less than a solar mass even in the absence of molecular cooling. This implies that low- 
mass primordial star formation is possible even without H2 if fragmentation of the clouds occurs 
at sufficiently high density. 

The outline of this paper is as follows. In §2, the method of our calculations is described. In 
§3, results of our calculations are presented. Finally, our work is summarized and its implications 
sketched briefly in §4. 



2. Model 

We consider spherical clouds of hydrogen-helium gas irradiated with FUV radiation. Specifi- 
cally, we investigate the cases of two FUV spectra in this paper: type a, the power-law radiation 

■AjvM = J21 x lO^VMh)" 1 ( er § s_1 cm ~ 2 str_1 Hz_1 ) (" < "th), (1) 
and type b, the diluted thermal radiation of 10 4 K 

Juv(^) = J21 x lO" 21 ^ 1 !!!^ (^g s" 1 cm" 2 str' 1 Hz' 1 ) [y < u th ), (2) 

where f t ii is the Lyman limit frequency. In both cases, JtjvC^) = beyond the Lyman limit. We 
neglect the so-called sawtooth modulation due to the absorption of Lyman series photons by the 
intergalactic matter (e.g., Haiman et al. 1997), for simplicity. This suffices for our simple analysis, 
which is valid only in order-of- magnitude estimates. We also include the cosmic microwave back- 
ground radiation (CMBR), although it has little influence on the matter at those low redshifts. The 
spectra of FUV radiation are shown in Figure 1 for J21 = 1 along with the CMBR at z = 30. Note 
that, for the same value of intensity at the Lyman limit, the number of photons above the threshold 
of H~ photodissociation is larger for type b spectra than for type a spectra. Consequently, the pho- 
todissociation rate of H~ is about 250 times larger for the T = 10 4 K thermal type radiation (type 
b) than for the power-law type (type a) with the same value of J21. Similarly, the photodissociation 
rate of H2 is about 3 times larger for the type b spectrum. 

We take the maximum expansion of an overdensity as the initial condition of the calculation. 
Assuming the redshift of that epoch z m ~ 30, we take the number density of hydrogen nuclei 
n = (f2bA^o)[p/(l + %He)wiH]) ionization degree y(e), and matter temperature T m at the maximum 
expansion epoch as n = 8.9 x 10~ 2 cm~ 3 ,T = 39K,y(e) = 2.0 x 10~ 4 , respectively. We neglect 
the initial molecular abundance and assume ^(Ehj) = 0. Here yne is the concentration of helium 
nuclei, and uih is the mass of a hydrogen nucleus. Note that the concentration of He is defined 
by UHe = nue/n, where n and nu e are the number density of hydrogen nuclei and helium nuclei, 
respectively. Simultaneously, we write for each atomic, molecular, or ionic species, y(x) = n(x)/n, 
where n(x) is the number density of species x. Note y(H.2) = 1/2 for fully molecular gas. We set the 
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primordial helium abundance at yn e = 0.0807, which corresponds to the mass fraction Y p = 0.244 
(Izotov k, Thuan 1998). The cosmological parameters are £Iq = l>^b = 0.05, and h = 0.7. 

Any effect due to rotation or magnetic fields are neglected for simplicity. Then, the actual 
collapse is expected to proceed like the Penston-Larson similarity solution (Penston 1969; Larson 
1969). 2 According to this solution, the cloud consists of two parts, that is, a central core region 
and an envelope. The central core region has a flat density distribution, whereas in the envelope, 
the density decreases outward as oc r~ 2 . The size of the central flat region is roughly given by the 
local Jeans length Aj = irc s / \J~Gp, where c s and p are the sound speed and total density in the core, 
respectively. In this paper, we take the radius of the central region as R c = Aj/2. The collapse 
in the core proceeds approximately at the free-fall rate. In this paper, we focus on the evolution 
in the central region. In particular, we calculate the the temperature and chemical composition of 
the collapsing core as a function of the central density. 

The cloud consists of two components, namely, baryon and dark matter. We assume that the 
dynamics of the baryonic component is described by the relation 

d Pb _ Pb 

dt {6) 
where p^ is the baryonic density in the central region and the free-fall time is 



Ws!?- (4) 

The dynamics of the dark matter is described by the relation for the top-hat overdensity, 

97T 2 , 1 + z„ 



m \3 



2 v l-cos6>' 



^oPc (5) 



up to the virialization, where the maximum expansion redshift z m = 30, and the parameter 9 is 
related to z by the relation 

l + . = (l + . m )(^^) 2/3 (6) 

7T 

(e.g., Padmanabhan 1993). We also use the usual age-redshift relation for the matter dominant 
universe, t = 3.1 x 10 17 (s)/i _1 f2 1//2 (1 + z)~ 3 / 2 . After the virialization of dark matter, i.e., when 
the density reaches to the virial density pvM,vir = 8pDM(zm), we keep pdm = PDM,vir- 

The thermal evolution is followed by solving the energy equation 

de d 1 A nct 

-rr = -P^;(— ) . ( 7 ) 

dt dt pb Pb 



2 Although the original Penston-Larson similarity solution is limited to the isothermal collapse, Yahil (1983) ex- 
tended this solution for general polytropic equations of state. 
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where 

1 kT 

e = (8) 

7ad - 1 |im H 

is the internal energy per unit mass of baryon, 

(9) 

limn 

is the pressure for an ideal gas, 7 a d is the adiabatic exponent, T is the temperature, \x is the mean 
molecular weight, ran is the mass of a hydrogen nucleus, and A nct is the net cooling rate per unit 
volume. The net cooling rate A nct consists of contributions from radiative cooling or heating by 
atomic hydrogen lines Ah; by rovibrational lines of molecular hydrogen, Ah 2 ; by continuum, A cont ; 
and by the Compton coupling with the radiation, Acjompt (the Compton cooling is unimportant 
in those low redshifts although it is included in our calculation) and from heating and cooling 
associated with chemical reactions, A cncm : 

Anct = A H + A Ha + A cont + Acompt + Achem- (10) 

The continuum processes included are listed in Table 1. The details of these processes are described 
in Appendix B. Since we are focusing on the evolution of the central region, whose radius R c = Aj/2, 
then we evaluate the optical depth t v by 

T v = K U R C = Kj,(y). (11) 

Time-dependent nonequilibrium chemical reactions are solved for the following nine species: 
H, H2, e, H + , H^", H , He, He + , and He ++ . Chemical reactions included are listed in Table 2. 



3. Results 

In this section, we present our numerical results. Figures 2 and 3 display the temperature 
evolution for collapsing primordial clouds irradiated with (1) the power-law type and (2) diluted 
black body of 10 4 K type FUV radiation, respectively. As mentioned in §2, the photodissociation 
rate coefficients of both and H2 are larger for type b spectrum than for the type a with the same 
value of J21. Consequently, FUV radiation of type b has the larger influence on the evolution of the 
clouds than type a with the same value of J21. We discuss effects of different spectral types and 
scaling relations between them in Appendix A. Except for those scalings of J21, the evolutionary 
trajectories change in similar ways with an increase of the UV intensity in both cases. Hence, we 
describe only the evolutionary features of the clouds irradiated with type a radiation in detail here. 

The evolutionary trajectories for the clouds with J21 > 10 5 are identical to each other. As is 
obvious from Figure 2, the thermal evolution at high densities (say, n > 10 7 cm~ 3 ) is completely 
different between clouds with J21 < 10 4 and those with J21 > 10 5 . This is because the clouds 
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irradiated by the FUV radiation with J21 > 10 5 cannot form sufficient H2, so they cool only by 
atomic cooling. The temperature of such clouds is higher than that of clouds with J21 < 10 4 , which 
can form enough H2 to cool eventually Toward higher densities, these two groups of trajectories in 
Figure 1 converge respectively to two different limiting tracks. We call the higher temperature one 
the "atomic cooling track" , and the other the "molecular cooling track" . This nomenclature comes 
from the fact that clouds evolving along the former track remain atomic, while those following the 
latter eventually become fully molecular. We should note that the Jeans masses inside the clouds 
ultimately fall below 0.1M Q for both tracks. 

Figure 4 shows the fractional abundances of hydrogen molecules (Fig.4a) and electrons (Fig.46) 
for the clouds irradiated with type a radiation. The separate heating and cooling rates per unit 
mass for those clouds are illustrated in the panels of Figure 5 for J21 = 0(Fig.5a), 10 4 (Fig.56) and 
10 5 (Fig.5c). 

We review their evolutionary features starting from the case of J21 = 0. In this case, the 
collapse proceeds along the well-known molecular cooling track (e.g., Palla et al. 1983). At the 
beginning, the temperature rises adiabatically owing to the compression and the initial lack of 
coolant. If the total mass of the cloud is less than the maximum Jeans mass attained at the end of 
this adiabatic phase, the contraction stops at the time when the cloud mass becomes equal to the 
Jeans mass. We assume here that the total mass of the cloud is larger than the maximum Jeans 
mass. At low densities (i.e., n < 10 8 cm -3 ), H2 is formed mainly through the channel: 

H + e -> H~ + 7 (reaction 32), (12) 

and then 

H + H~ -> H 2 + e (reaction 9). (13) 

Increased temperature and density causes an increase of both the H2 formation rate and the H2 
cooling function. Consequently, a sufficient amount of H2 to cool within a free-fall time is formed, 
and the temperature then drops. A detailed discussion of when the efficient H2 cooling begins is 
presented in Appendix A. Owing to the H2 cooling, the temperature is kept as low as a few hundred 
K in the cloud (see Fig.2 and Fig.5a). When the density reaches about 10 8 cm~ 3 , the three-body 
reaction 

3H -> H 2 + H (reaction 16) (14) 

becomes efficient (Palla et al. 1983). In consequence of this highly productive H2 formation, 
hydrogen becomes fully molecular (see Fig.4 a). At the same time, the cloud becomes optically 
thick to a few H2 lines at n ~ 10 n cm~ 3 . However, the cooling remains efficient enough to induce 
dynamical collapse (i.e., The ratio of specific heat V = dlogp/dlogp < 4/3) because there are still 
enough optically thin lines. When the central number density reaches about 10 14 cm -3 , H2 collision- 
induced emission begins to dominate the cooling (see Fig.5a). The cloud becomes optically thick to 
this continuum at ~ 10 16 cm~ 3 . Thereafter, the radiative cooling rate drops rapidly. Simultaneously, 
H-2 dissociation begins, because the temperature is already near the value needed for dissociation. 
The H-2 dissociation prevents the temperature from rising rapidly until the number density reaches 
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about 10 20 cm~ 3 . The minimum Jeans mass of 1.5 x 10~ 3 M© is reached shortly after that. When 
most of the hydrogen molecules are dissociated, the collapse becomes approximately adiabatic. 
According to Omukai & Nishi (1998), after a little further contraction, a small hydrostatic core 
(about a few 1O -3 M0), or protostar, forms at n ~ 10 22 cm~ 3 ,T ~ 3 x 10 4 K. 

As seen in Figure 2, if the intensity of FUV radiation increases, the density and temperature 
at the time when efficient H2 cooling begins also increases (see Appendix A for further discussion). 
After H2 cooling becomes effective, the trajectories of those clouds soon converge to the evolutionary 
track for J\jy(v) = 0. 

An interesting example is when the intensity is as high as J21 — 10 4 . In this case, the temper- 
ature rises almost adiabatically up to about 8000K without enough H2 formation. At that point, 
Lya line emission of atomic hydrogen begins to work as an efficient cooling agent (see Fig. 5 b). After 
the cloud collapses isothermally at T ~ 8000K over an order of magnitude in density owing to this 
coolant, sufficient hydrogen molecules form as a result of the increased density. Thereafter, the 
trajectory converges rapidly to the molecular cooling track in the same way as those with weaker 
radiation (see Fig.2). 

When the intensity of FUV radiation is increased more, i.e., J21 > 10 5 , the thermal evolution 
becomes completely different from those with lower intensities (see Fig.2). The adiabatic rise of 
temperature continues until Lya cooling becomes effective at about 8000K. 3 While the cloud 
collapses nearly isothermally at about 8000K owing to Lya line cooling, the central density reaches 
the critical density 4 n cr ~ 10 4 cm~ 3 of hydrogen molecules before a amount of H2 sufficient to cool 
within a free-fall time forms. At higher densities, n > n cr , enough H2 never forms, for the following 
two reasons. First, the H2 fraction decreases; the collisional dissociation rate coefficient increases 
because, at n > n cr , high vibrational levels, at which dissociation occurs easily, are populated by 
more hydrogen molecules than at lower density. Second, the amount of H2 needed for efficient 
cooling increases; note that the cooling rate per unit volume Ah 2 oc n 2 for n < n cr , while Ah 2 oc n 
for n > n cr . (see Appendix A for further discussion). Consequently, the cloud never joins to 
the molecular cooling track. Instead, it continues to cool by Lya line emission until the density 
reaches about 10 6 cm -3 , where collisional de-excitation from the 2p state begins to dominate the 
Lya emission as a result of the small escape probability (i.e., /3 esc ,2i^42i — C21). The subsequent 
major cooling process is two-photon emission of atomic hydrogen from the 2s state. This works 
effectively in the density range 10 6 cm -3 < n < 10 7 cm~ 3 . Next, the cloud cools by the free-bound 
emission of H~ over nearly ten orders of magnitude in density, until the number density reaches 
10 16 cm -3 (see Fig. 5 c). The mechanism of H free-bound emission cooling is as follows. As a first 



3 The initial slight rise of ionization degree seen in Figure 46 is caused by photoionization from the first excited level, 
which is populated by more atoms than the local thermodynamic equilibrium (LTE) value as a result of absorption 
of irradiated Lya photons. This effect has little influence on the later evolution. 

4 The critical density is the density above which the rovibrational levels of H 2 are populated, according to the LTE 
law. 
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step, radiative association of H and e; 

H + e -> H _ + 7 (reaction 32) (15) 

occurs and a photon is emitted. This H is used to form H2 (reaction 9), which will be collisionally 
dissociated by reaction 15. As a result, there is net cooing by the photon emitted in reaction 32. 
5 In the course of the dynamical collapse induced by H~ free-bound emission, the temperature 
falls gradually from 7000K (at n ~ lOW" 3 ) to 3000K (at n ~ 10 16 cm- 3 ; see Fig. 2). At the 
central density, ~ 10 16 cm~ 3 , the cloud becomes optically thick to both H~ bound-free absorption 
and Rayleigh scattering of atomic hydrogen. At the number density ~ 10 17 cm -3 , the H ionization 
begins (see Fig. 4 b). As a result, the temperature rise slows down up to n ~ 10 19 cm~~ 3 . At 
n = 10 20 cm , where 60 % of the gas is ionized, the minimum Jeans mass Mj iin j n = 0.03Mq is 
reached. Thereafter, the temperature continues to rise almost adiabatically. We expect that a 
protostar on the order of Mj >m - m is formed after a little further adiabatic contraction. Although the 
minimum Jeans mass of the atomic cooling clouds is 20 times larger than that of the H 2 cooling 
clouds, it is still much less than the mass of usual stars. 



4. Summary and Discussion 

We have investigated the thermal and chemical evolution of primordial clouds irradiated with 
FUV radiation. When intensity of the irradiating FUV radiation exceeds some critical value [J21 — 
10 5 for Jtjv(^) °t 1/1 (hv < 13.6eV)], sufficient molecular hydrogen to be important for cooling is 
never formed in those clouds because of the photodissociation and blocking of the formation channel. 
Nonetheless, sufficiently massive clouds can start dynamical collapse by atomic line cooling even 
without H2. Those clouds continue to collapse almost isothermally at several thousand K along 
the "atomic cooling track" in the density-temperature plane. That dynamical collapse is induced 
by successive cooling by Lya emission, two-photon emission, and H free-bound emission. The 
minimum Jeans mass eventually reduces to about O.O3M0 for those clouds. On the other hand, 
clouds irradiated with less intense FUV radiation collapse dynamically by the H2 cooling. The 
minimum Jeans mass of such clouds is about 1.5 x 1O~ 3 M . 

Although in this paper we limited our application to the protostellar collapse under an FUV 
background, we expect that the primordial protostellar collapse without H2 also occurs in other 
situations. Recently, Susa & Kitayama (2000) investigated the collapse of primordial clouds ir- 
radiated with ionizing UV background Jtjv(^) = J21 x 1® 1 {v / vth) -1 (erg s _1 cm^ 2 str _1 Hz _1 ). 
They pointed out that H2 formation is suppressed for J21 > 3.6 x 10 4 owing to blocking of the 
H2 formation channel by H~ photodissociation. Although they did not study the later evolution, 
we guess that the clouds continue collapsing along our atomic cooling track. Carlberg (1981) also 



5 Formed H is photodissociated at low densities, that is, n < k^/kg = 0.1a J21, where a is a constant defined in 
Appendix A, and is 8 for the type a (2 x 10 3 for type b, respectively) spectrum. In this case, there is no net cooling. 
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considered the thermal evolution of collapsing primordial clouds without external radiation, but 
he did not include the three-body process of H2 formation. For this reason, in his calculation, 
molecular hydrogen dissociates at n ~ 10 n cm~ 3 . After that, the evolutionary path of his cloud 
follows approximately our atomic cooling track (see Fig. 2 of Carlberg 1981). This fact support 
our speculation that the clouds that do not cool by H 2 collapse along the atomic cooling track. 

Next, we discuss the feasibility of the collapse along the atomic cooling track in the cosmological 
reionization. The UV intensity needed to reionize the universe at z re \ OD , i.e., 1 ionizing photon per 
a hydrogen nucleus, is 

Jrcion = k f b ,f C f A + Z : ei ° n)3 = 3 x l(n 21 (i±^f (erg s- 1 cm- 2 str" 1 Hz" 1 ) (16) 
47r(l + 4y Hc )"iH 10 

at the Lyman limit for J{y) oc v~ l and the assumed cosmological parameters. Thus, the condition 
that the clouds collapse along the atomic cooling track, J21 > 10 5 for Juv(^) °^ h seems to be 
rarely met in practice, except in the immediate vicinity of UV sources. 

Finally, we discuss the mass of formed stars. Although the minimum Jeans mass is less than 
0.1 Mq, this corresponds not to the final mass of formed stars, but rather to the initial mass 
of the protostar, which has no direct relation with the former. When a small hydrostatic core, 
namely protostar, forms at the center of a collapsing protostellar cloud, an enormous amount 
of gravitationally unstable gas still surrounds the protostar. Thereafter, the protostar grows in 
mass by accreting the ambient matter. Consequently, the final mass of stars is determined by the 
subsequent accretion onto the protostar. Stahler, Shu, & Taam (1980) argued that a rough estimate 
of the protostellar mass accretion rate M can be obtained from the relation M ~ c 3 /G, where c s 
is the isothermal sound speed in the protostellar cloud. From this relation, the mass accretion 
rate is higher for protostars formed by atomic cooling than for those formed by molecular cooling 
because of the higher temperature of the protostellar cloud. Suppose that the accretion is halted 
by some stellar activity, e.g., radiation force, bipolar flow, etc. In this case, the higher accretion 
rate probably results in the later halting and higher mass of formed stars (e.g., Wolfire & Cassinelli 
1987). 

In contrast, from the viewpoint of fragmentation, the masses of primordial stars formed by 
atomic cooling alone could be smaller than those formed by molecular cooling. The fragmentation 
is possible in principle as long as the temperature decreases toward higher density, while it does not 
occur after the isothermality of the cloud breaks down and the temperature begins to rise (Inutsuka 
& Miyama 1997; Masunaga & Inutsuka 1999). From Figures 2 and 3, the final fragmentation occurs 
at about n ~ 10 16 cm -3 for clouds collapsing along the atomic cooling track. On the other hand, it 
seems difficult for molecular cooling clouds to fragment at n > 10 3 cm -3 . The masses of fragments, 
which is about the Jeans mass at that time, are about O.1M for the atomic-cooling clouds, and 
1O 4 M for the molecular-cooling clouds. If this is the case and the whole material contained in a 
fragment is converted into a single star, very massive stars form from a molecular-cooling cloud, and 
red dwarfs form as a result of atomic cooling. However, it should be noted that the fragmentation 
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and thermal evolution depend strongly on each other and also on the geometry of the clouds. In 
fact, Uehara et al. (1996) concluded that the minimum mass of fragments is on the order of a solar 
mass, which is essentially the Chandrasekhar mass, by studying thermal evolution of filamentary 
primordial clouds that collapse by molecular cooling. Later, improved simulations by Nakamura 
& Umemura (1998) confirmed Uehara et al.(1996)'s result. The discrepancy between our mass of 
fragments (~ 10 4 M Q ) and theirs (~ 1M Q ) results from the difference of the assumed geometry of 
the clouds; we have assumed spherical clouds in studying the thermal evolution, while their clouds 
are filamentary. For filamentary clouds, the gravitational contraction is slower and compressional 
heating rate is lower than for spherical ones. Hence, for the filamentary cloud the isothermality 
breaks down at higher density, where the Jeans mass is smaller. To fully address the complexity 
arising from the cloud geometry, two-dimensional or three-dimensional calculations are needed. 
Furthermore, there are other uncertainties relating to the evolution of the fragments, for example, 
mass accretion onto the fragments, merger between them, feedback from neighboring formed stars, 
etc. (e.g., Bromm et al. 1999). It seems premature to draw conclusions here about these issues. 

I would like to thank Shu-ichiro Inutsuka, Ryoichi Nishi, Humitaka Sato, Toru Tsuribe for 
crucial suggestions and discussions, and Phil Stancil for providing the data of the H^" photodisso- 
ciation cross section. I also acknowledge the referee for improving the manuscript. This work is 
supported in part by Research Fellowships of the Japan Society for the Promotion of Science for 
Young Scientists, grant 4287. 
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A. Conditions for H2 Cooling 

In this section, we discuss key processes determining the amount of the molecular hydrogen 
forming under FUV radiation. In particular, we delineate the condition for sufficient H2 for cooling 
to form. 

The H2 fraction i/(H 2 ) necessary to cool within a free-fall time is given by 

(3/2)fc B T 

2/cooi(H 2 ) = — - — — , (Al) 
nLB_ 2 ts 

where Ln 2 (erg s -1 cm 3 ) = AH 2 /n(H)n(H2) is the H2 cooling function and we assume n(H) ~ n. 
When the actual H 2 fraction y(H2) reaches y C ooi(H2), the temperature drops abruptly as seen in 
Figure 2 and 3. 

Next, we discuss how much H2 forms during the collapse. The H2 abundance is bounded by 
the following three values; 

(1) the amount of H2 that is able to form within a free-fall time, yf or m(H2), 

(2) the chemical equilibrium value between H2 formation and photodissociation, y p d(H2), 

(3) the chemical equilibrium value between H2 formation and collisional dissociation, y c d(H2). 
In practice, the smallest among these three values is realized; 

y(H 2 ) =min[2/ form (H 2 ),yp d (H 2 ),ycd(H2)]. (A2) 

In the following, we examine the values of yform(H2), y p d(H2), and y c d(H2) in this order. 

First, we find the formation time limited value yform(H2). H2 formation occurs mainly through 
the H~ channel: 

H + e -> H~ + 7 (Reaction 32), (A3) 

followed by 

H + FT -> H 2 + e (Reaction 9). (A4) 

In the above chain, not all of H~ formed are used in H2 formation. Instead, some of H~ arc 
photodissociated via Reaction 31; 

FT + 7 -► H + e (Reaction 31). (A5) 

The photodissociation rate coefficient of H - , ^31, depends on the spectrum of radiation. Here, we 
introduce a parameter a, defined by 

hi = (XK31J21, (A6) 

where K31 = 1 x 10 -10 is the value of k%i under Jtjv(^) = 10 -21 = const, (for v < u t ^). The 
parameter a characterizing the radiation spectrum above H~ photodissociation threshold (0.755 
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eV) is 8 for the power-law type spectrum (type a in the text), 2 x 10 3 for the diluted thermal type 
spectrum (type b). The rates of the competing reactions 9 and 31 determine the branching ratio 
of H _ to be used in H2 formation, which is kgn/(kgn + &3i). Then, by defining the effective H2 
formation rate coefficient 

fcform = h2-, —T~ > ( A7 ) 

Kgn + K31 

the H2 formation rate can be written as fcf orm n(e)n(H). Using kf orra above, we obtain the H2 fraction 
that can form within a free-fall time 

yform(H 2 ) = k {oTia y(e)nt s . (A8) 

Second, we find the photodissociation limited value y pd (H 2 ). The photodissociation occurs via 

H 2 + 7 — ► H| — ► 2H (Reaction 35), (A9) 

whose reaction coefficient is 

k pd = k 35 = ltfj(hv = 12.4eV)/ sh = 10- 12 /3J 2 i/ sh . (A10) 

Here, we have introduced a parameter (3 that represents the relative strength of radiation intensity 
at the average LW band frequency (12.4 eV) to that at the Lyman limit (13.6 eV). The shielding 
factor / s h is given by (Draine & Bertoldi 1996) 

f sh = min [l, (AWl0 14 cm- 2 )- 3 / 4 ] • (All) 
The equilibrium value between the H2 formation and photodissociation is 

y P d(H 2 ) = %^y(e)n. (A12) 

If Nu 2 > 10 14 cm~ 2 , the coefficient &pd, and then the right hand side of equation (A12), involves 
y(H 2 ) itself. Denoting k pd = a pd y(H 2 )~ 3/4 , where a pd = 7 x 10 _18 T 3 ~ 3/8 n~ 3 / 8 /3J2i, we obtain 

y P d(H 2 ) = (^y(e)n) 4 (N H . 2 > 10 14 c m - 2 ). (A13) 

Combining expressions for cases with and without H2 self-shielding, the photodissociation limited 
value of H2 is 



y P d(H 2 ) = max 



fcform 1 \ ( ^form 



y(e)n, y(e)n 



(A14) 



10- 12 /3J 2 i v ; ' V « P d 

Note that the radiation intensity enters into the expression for y pd (H 2 ) only in the form of afiJ^i 
in the low density range n < k^i/kg ~ O.la^icm" 
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Third, the collisional dissociation limited value y c d(H2) can be obtained as follows. The dom- 
inant collisional dissociation process is collision with the atomic hydrogen: 

H 2 + H -> 3H (Reaction 15). (A15) 
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Then, denoting /c C( j = k±5, 

y cd (H 2 ) = ^y(c). (Ai6) 

This value y c d(H2) becomes smaller when the density exceeds the critical density n CT ~ 10 4 cm~ 3 
owing to the collisional dissociation from high vibrational states. 

Now, we examine when sufficient H2 forms under FUV radiation in the case that the tem- 
perature rises adiabatically from the initial value. Figure 6 shows the H2 fraction estimated by 
Equation (A2) in this case. Here, the temperature is given by T = To(n/no) 2//3 , where the initial 
number density and temperature are no = 8.9 x 10~ 2 cm _3 ,To = 39K, and the ionization degree 
y(e) is fixed to 10" 4 . 6 In Figure 6, each curve for y(H2) with fixed FUV intensity consists of four 
segments. From the lowest to highest temperature in this Figure, those segments correspond to 
y p d(H2) without self-shielding (gradually increasing portion), y p d(H2) with self-shielding (rapidly 
increasing portion), yf or m(H2) (another gradually increasing portion), and y c d(H2) (decreasing por- 
tion). As known from Figure 6, except for lowest intensities (in such a case, the evolution of the 
cloud is virtually unaffected by the presence of FUV radiation) , the photodissociation limited value 
ypd(H2) determines when the sufficient H2 forms. The condition that molecular cooling becomes 
effective before atomic cooling does is that y p d(H2) > y CO oi(H2) at T = 8000 K, which leads to 

J 2 i<4xl0 3 (^) 1 /2, (A17) 

where X4 = y(e)/10 -4 and we have used kf orm ~ k^kgn/k^i, and /C32 = 3.6 x 10~ 15 ,/cg = 8.7 x 
10~ 10 ,y C ooi(H2) = 10~ 6 , and n = 2.6 x 10 2 cm" 3 at 8000K. Also, since FUV intensity enters into 
y p d(H2) in the form of cc/SJfi, type b spectrum of has the same effect as type a spectrum of 
i/a/3(b)/a/3(a) ~ 27 times on the early evolution of the clouds. 

As discusses in the text, the cloud irradiated by strong FUV radiation never forms sufficient 
H2 for cooling if the intensity exceeds a certain critical value. Next, we discuss the physical reason 
where this value comes from. Here, we examine how much H2 forms while the cloud collapses 
isothermally owing to Lya cooling. The values of H2 concentration y(H 2 ) estimated by equation 
(A2) for T = 8000K and y(e) = 10~ 4 are shown in Figure 7. 7 For a given value of FUV intensity, 
y(H2) is limited by y p d(H2) (an almost vertical line) in the lower density range and by y c d(H2) (a 
smooth curve) in the higher density range. The formation time limited value J/f rm(H2) does not 
appear in this Figure. On the other hand, y C ooi(H2) decreases in the low density regime, reaches 
the minimum at n ~ n CI ~ 10 4 cm~ 3 , and increase in higher densities. As seen in Figure 7, the 
H2 fraction y(H 2 ) has no chance to reach y coo i(H2) in n > n cr — 10 4 cm -3 . This is because the H2 



Here and through out this section, we assume the ionization degree is constant, namely, the recombination 
and/or ionization time scale is longer than the evolutionary timescale of the system (i.e., free-fall time) Although this 
is somewhat oversimplification, it suffices for clarifying key processes determining the H2 fraction. 

7 Once the actual H2 fraction exceeds t/ C ooi(H2), the temperature of the cloud drops from ~ 8000K as a result of 
the H2 cooling. After that, the H2 fraction shown in the Figure does not represent the actual value. 
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fraction y(H-2) is limited by y c d(H2), which becomes small in those densities owing to the enhanced 
collisional dissociation from high vibrational levels. Therefore, in order for the H2 cooling to become 
efficient, y(H.2) must reach y C ooi(H2) at n < n cr ; both y c d(H2) and y p d(H2) must exceed y C ooi(H2) 
at n < n cr ~ 10 4 cm -3 . As the condition on y c d(H2), we take "fy c d > y c d(J(v) = 0), where constant 
7 ~ several, noting that y c d(H2; J{y) = 0) is several times larger than y C ooi(H2) at n < n CT . This 
reduces to 7fcgn cr > k^i, or in another expression, 

J21 < 4 x lr/cT 1 , (A18) 

where we have chosen 7 = 5. The condition on y p d(H2) leads to another constraint on J21: 

J21 < 1 x W 6 x 4 (3~ l . (A19) 

Here, we have used y coo i — 1CT 6 , since the dependence of y p d(H2) on n is so strong that the precise 
value of y C ooi(H2) is not important. Also, fcf orm = k%2 has been used. Combining the conditions 
(A18) and (A19), the critical intensity J cr above which the cloud never forms sufficient H2 is 

J cr 21 = min(4 x 10 5 a-\ 1 x 10 6 x 4/ 9 _1 ). (A20) 

Substituting the parameters a and (3 into equation ( A20) , we see that the former constraint is more 
restrictive than the latter for both two types of radiation discussed in the text. 



B. Radiative Processes 

B.l. Atomic Hydrogen Lines and Two-photon Emission 

We model an atomic hydrogen as a five level system. The radiative cooling rate due to atomic 
hydrogen lines per unit volume is given by 

Ah lines = ^2 hv u i(3 cscM A u in u (B)[SB.{v u l) - J C ont{^ul)]/Su{^ul), (Bl) 
ul 

where n u (H) is the population density of atomic hydrogen in the upper energy level u, A u \ is the 
Einstein radiation probability for a spontaneous downward transition, /3 CS c,n/ is the probability for 
a emitted line photon to escape without absorption, and hv u i is the energy difference between the 
upper level u and the lower level Z, and Jcoat{ v vl) is the mean intensity of overlapping continuum 
at the line frequency v u \. The source function Su(v u i) is given by 

5hM "— ( ^H) } • (B2) 
The cooling rate due to the two photon emission is 

A 2p h = /i^2i/3csc,2 P hA 2 in2(H), (B3) 



-15- 



where A21 is the spontaneous downward transition probability by the two-photon emission. The 
escape probability /3 C sc,2ph = ex P( — r a)> where r a is the absorption optical depth owing to other 
continuum processes (al-a7 in Table 1; See §B.3). We take only account of absorption by other 
continuum processes, since our clouds do not become optically thick to the two-photon continuum. 

Relative population of each energy level is obtained from a solution of the equations of the 
detailed balance between levels (e.g., de Jong, Dalgarno, h Boland 1980; Tielens & Hollenbach 
1985) 

n n 

m(B) Rij = TuQVRji, (B4) 

where n is the total number of transitions included. The transition rate Rij from level % to level j 
is given by 

ft.. = f {AijPescij + A ijA;sc,2ph)(l + Qij) + Qj for 1 > j 

\ (9j/9i)(^-ji(3csc,ij + ^jil3 CS c,2ph)Qji + Cij for % < j, 
where Cij is the collisional transition rate, and 

c 2 

Qij = 2hv % Jcmit{nj)- (B6) 

The relative population within the first excited states (i.e.,2p and 2s states) is obtained from 
the statistical equilibrium (Spitzer & Greenstein 1951), 

n 2s _ 92s ^ C 2s 2p ^ j-g^ 
n 2p 92p Cis2p + ^2sls 

where the statistical weight g2 S = 2, and g2 P = 6, the radiative transition rate by two photon 
emission A2 S i s = 8.23(s~ 1 ), and the collisional transition rate between the levels 

C 2s2p = 6.21 x 10- 4 T- 1 /2] n (5.7T)[l + J -^ ?ff - ) }n(e)(s~ 1 ). (B8) 



Using ri2p obtained above, we may write 



A 21 = r ^A 2 p ls (B9) 
™2 



where A 2p i s = 6.27 x lO 8 ^ 1 ). 



We assume LTE within levels of the same principal quantum number for n > 3. Averaged 
over angular momentum quantum numbers, A31 = 5.575 x lO 7 ,^! = 1.278 x 10 7 ,yl5i = 4.125 x 
10 6 ,A 32 = 4.410 x 10 7 ,^42 = 8.419 x 10 6 ,A 52 = 2.530 x 10 6 ,^ 4 3 = 8.986 x 10 6 ,A 53 = 2.201 x 
10 6 , A 54 = 2.699 x 10 6 (Janev et al. 1987). 

Likewise, we use the two-photon emission rate A21 = {n 2s /n 2 )A2 S i s and A n ; = for other 
transitions. 
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The collisional de-excitation rate is given by 

Cui = 7 „/(e)n(e) + 7lli (H)n(H), (BIO) 

where the collisional de-excitation rate coefficients 7 u ;(e) for collisions with electron and j u i(H) for 
those with atomic hydrogen are given below. 

The collisional de-excitation rate coefficients for collisions with electron is given by 



= 10-(^)3/4 ai „i|2±S, _ (B11) 

(Sobelman et al. 1981), where ai 2 = 24, a 13 = 22,a 14 = 22,a 15 = 21,a 23 = 67, a 2 4 = 58, 025 = 
56,a 34 = 124, a 35 = 101, 045 = 185 and X n = 0.28, X i3 = 0.37, X i4 = 0.39, X is = 0.41, X 23 = 
0.30, X 24 = 0.45, X 25 = 0.52, X34 = 0.26, X 35 = 0.42, X 45 = 0.21. For collisions with atomic hydrogen 
(Drawin 1969), 

, s „ 9 , 9x 9 „ 1/9 1 + 1.27 x 10~ 5 (1// 2 - l/n 2 )-^ 

7„i(H) = 7.86 x 10" 15 (- 2 (1// 2 - 1/u ) fi u T / — — ' ' „ „ , (B12 

' V 7 V v/ 77 1 + 4.76 x 10- 17 (1// 2 - l/n 2 )" 2 T 2 ' v 7 

where / i2 = 0.4162, f 13 = 7.910 x 10" 2 , f 14 = 2.899 x 10~ 2 , f 15 = 1.394 x 10" 2 , / 23 = 0.64 07, / 24 = 
0.1193, / 25 = 4.467 x 10~ 2 ,/ 3 4 = 0.8421, / 35 = 0.1506, / 45 = 1.038. 

The collisional excitation rate can be obtained from the detailed balance: 

Cm = C ul {g u / gi )exp(-hu u i/kT). (B13) 

Taking into account the collisional de-excitation of line photons and absorption due to the 
overlapping continuum, we may write the escape probability as 

rn^' CSC 

Pesc,ul = 2' exp(-T a ), (B14) 

where p u i is the probability for an absorbed line photon to re-emerge as an original line photon 
with neither collisional de-excitation nor two photon continuum emission, and N csc is the number 
of scattering that an average photon experiences before escape. The re-emergence probability can 
be written as 

p ^ui (B15) 

A u i + C u i + A u i 

The average number of scattering for Lya line in large optical depth r 2 i limit is given by (Bonilha 
et al. 1979) 

A^csc = 1.675i(l + O.OSC 1 - 5 )- 1 (B16) 

for escape from a moving infinite slab, where £ = 0.90t>b u ik/^D, ^buik is the velocity difference 
between the center and the edge of the slab, and = ^2kT /m-R. Here, as Wbulk, we take the 
velocity at the edge of the central core region whose radius is R c . Then Wbuik = Rc/^tff- 
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For other atomic lines, we take 

f l.Gr ul forr ui <10 4 
[ 160rJ for t u[ > 10 4 

(e.g., Elitzur & Ferland 1986). The optical depth averaged over a line is given by 

A u ic 3 

T u i = -^riMtygu/gi - n u (R)]l sh /vT). (B18) 

In order to take into account the effect of velocity gradient, we reduce the shielding length Z s h for 
lines other than Lya in the case of large velocity gradient as 

l sh = mm(R c , As th ) (B19) 

where 

civ 

Asth = W(^) = 3w D%. (B20) 

Recall that the correction for velocity gradient has been already included in the expression (B16) 
for Lya photons. Then we take 4h = Rc for the Lya line. 



B.2. Molecular Hydrogen Lines 

The radiative cooling rate Ah 2 due to molecular hydrogen lines can be represented by the 
similar expression as equation (Bl). We compute rovibrational population n t) j(H2) of H2 following 
Hollenbach & McKee (1979) using the renewed collisional rate coefficients given by Hollenbach &; 
McKee (1989). We take the spontaneous radiative decay rates A v jyji from Turner, Kirby-Docken, 
& Dalgarno (1977). Also we take into account the effects of the external radiation by the same way 
as described for the atomic hydrogen lines (i.e., equation B5). The escape probability is given by 
(Takahashi, Hollenbach, k Silk 1983) 

/3csc,«/ = [(l-e-^)/^]e- Ta , (B21) 

for the case that the velocity is proportional to the radius. We consider the first three vibrational 
states (v = — 2) with rotational levels up to J = 20 in each vibrational state. We assume the 
ortho to para ratio of molecular hydrogen to be the equilibrium value 3:1. We use the rovibrational 
energies E(v, J) of Borysow, Frommhold, & Moraldi (1989). 



B.3. Compton Coupling 



The cooling rate of a gas with electron density n e and temperature T embedded in a black body 
radiation field of energy density u and temperature T 7 is given by (Kompaneets 1957; Weymann 
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1965) 

Acompt = — c2 7 caTUeU (ergs s^^cnT 3 ) (B22) 

= ^^ T^T-T,) (ergs s^cm" 3 ), (B23) 

where ut is the Thomson cross section, <tb is the Stefan-Boltzmann constant, and m e is the mass 
of an electron. In order to treat the situation that the Compton energy transfer and continuum 
processes coexist both in optically thin and thick cases, we deal with the Compton process as if it 
was one of continuum sources, denoting 

KCompt(^) = ^^0"T«e = OCompt^, (B24) 

and 

AkT 

f?Com P tO) = 7,a T n e J(u) = acomptTJ(v), (B25) 

m e c z 

where acompt = (4/c/m e c 2 )a"Tra e . This method reproduces the correct behavior both in optically 
thin and thick limit (Note that T 7 = T in the optically thick limit). In evaluating the radiation 
temperature in the above equation, we use the definition 

T 7 = (7rJ/a B ) 1/4 , (B26) 

where the mean intensity J = J J{y)dv. 



B.4. Continuum 

We consider the following continuum processes (cf., Lenzuni et al. 1991). Among pure ab- 
sorption, bound-free absorption of H, He, H~, H^", free-free absorption of H~, H, H 2 collision-induced 
absorption, 8 and among scattering, Rayleigh scattering of H, and Thomson scattering are included. 
The cross sections are listed in Table 1. 

The net rate of energy transport from matter to radiation per unit volume per unit frequency 

is 

AcontW = 47r[77» - K a (u)J(u)}, (B27) 

where ?/ a (z / ) is the thermal part of the emission coefficient and K a (^) is the true absorption coeffi- 
cient. Also, we use the scattering coefficient k s (v) below. Specifically, 

a8 

1aM= £ (B28) 

aj=al 



Two-photon emission of atomic hydrogen was treated as described in §A.l. 
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and 

a8 

K aH = ^ K a j(», K s (u) = K sl (u) +K s2 (v). (B29) 
aj=al 

We assume that the radiation field is static inside the cloud, then 

Acont(,) = (B30) 

where u(u) and u cx (u) are the radiation energy density in side and outside the cloud. The photon 
escape time 

wH = 4t-> ( B31 ) 

where the extinction coefficient 

XW = «aW+«sM, (B32) 

and 

x{v) = max[r(z/) 2 , t(is)], (B33) 

by using optical depth t(v) = x{ u )Rc- Equating equations (B27) and (B30), and using u{v) = 
AttJ(u)/c, we obtain 

J(u)= (B34) 

where = re a (zv)/x(z/). The external radiation field J cx (^) = •/uv( I/ ) + -B(i/;Tcmbr) i n our 
case, where Tcmbr is the temperature of the CMBR. If we include the Compton energy transfer, 
substituting 

?7aH +?7Compt(^) = Vaiy) + ClComptT J (l>) 
K a (l/) + KCompt(^) Ka(f) + OCompt^ 

into equation (B34) instead of S a (V), we obtain the final expression for the mean intensity inside 
the cloud: 

j, x = ^cx(^) + £(^)s( I/ ) T 7a(f)/(«a(f) +QComptT 7 ) 

1 + £(i/)x(i/)[l - acom P tT/( Ka H + a C omptT 7 )] ' 

Substituting equation (B36) into (B27) and integrating over frequency u, we obtain the sum 
of cooling rates A cont + Acjompt due to the continuum and the Compton coupling. 

In the following, we describe how to calculate the emission and absorption coefficients using 
the cross sections. 



B.4-1- bound-free absorption / free-bound emission (al-a4) 

Consider the radiative association RA(i) of species A and B into ith state of C, whose binding 
energy is hi/i, and resultant free-bound emission, 

RA(i) : A + B -> C(i) + hi/, (B37) 
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and its inverse reaction PD(«)(i.e., photodissociation of C from ith state, or the bound-free absorp- 
tion). The cross section of this (and its inverse) reaction is <7 R xi{v){a-p-oi{u) , respectively). In our 
case, the species C means H, He, H~, and for the processes al-a4, respectively. 

From the Milne relation (e.g., Rybicki & Lightman 1979), 

hv 2 2g Ci , , 

ORAi = crpDii ) , (B38) 

m T cv zazb 

where the reduced mass m T = m^m-Q/mQ, ^a(^b) is the partition function of A (B, respectively), 
gci is the statistical weight of the ith state of C, and 

^m r v 2 = h{y - Vi). (B39) 

The emission coefficient of RA(i) is given by 

riBAiiy) = ^a RAl n(A)n(B)vf(v)^ (B40) 

where the distribution function of the relative velocity between particles A and B 

Also, the absorption coefficient of PD(i) is given by 

KPDi(^) = cr PDi ni(C) - ^r^3°'RAjn(A)n(B)u/(t;)— (B43) 

= ™( c »-S ( ra )3/2w( -T 1 »"( A) "( B »- < B44 > 

The second term in the equation above represents the induced association. 

Summing over all levels, we obtain the emission and absorption coefficient owing to the reaction 

RA : A + B -> C + hv, (B45) 

and its inverse PD; 



2hl/ \* „tn\^/ ^u^W n (A)n( B 

17F ~^c) 



2 a PD n(C)exp(- — )[K(T)^ \>^\ ' ], (B47) 



and 



^pd(^) = J^KiW (B48) 

i 

= . PD „(C){l-e x p(-^, l ^A'(T)-!^W B l]}. (B49) 
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In the above equations, the equilibrium constant 

K(T) s ,»(%(a r (B50) 

_ 25,^,^,^,, (B51) 

and 

0"PD = ^2 CiCTpDi ' °"pd = ^2 c *°"PDi, (B52) 

i i 

where the relative population of level i is Cj = nj(C)/n(C), and its LTE value c* = gciexp(— — 
Vi)/kT)/z c . 



B.4-2. free-free absorption/ emission (a5,a6) and H2 collision-induced absorption/ emission 

(a7,a8) 

First, we consider the free-free emission and absorption (a5,a6) 

FF:A+e^A+e + 7, (B53) 

where A=H, and H + for the processes a5, and a6, respectively. The absorption cross section for A 
is off- Note that off oc n(e) (see Table 1). This process is collisional, then it occurs at the LTE 
rate. Thus, 

VFFiv) = —^-o- FF n(A)exp(-—), (B54) 

and 

kff(^) = a FF n(A){l - exp(-^)}. (B55) 

Next, we consider the H2 collision-induced emission and absorption (a7,a8) 

CI : H 2 + B ~ H 2 + B + 7, (B56) 

In this equation, B means H 2 for the process a7, and He for a8. The absorption cross section for 
H 2 is (Jci- This process occurs only at so high density that aci = Vqi- Then, in the same way as 
the above, 

r/Ci(f) = -^ 2 - cr ci«(H 2 )exp(-— ), (B57) 

and 

kciM = ^Cin(H 2 ){l - exp(-^)}. (B58) 
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B.4-3. scattering (sl,s2) 

We denote the cross section of A for the scattering (sl,s2) 

SC : A + 7 A + y (B59) 

as <rsc(^)- In the above expression, A=H, and e for the scattering si, and s2, respectively. The 
scattering coefficient is given by 

ksc(v) = <tsc(vMA). (B60) 

C. Chemical Reactions 

We treat non-equilibrium chemistry of hydrogen-helium gas between the following nine species: 
H, H2, e, H + , H^, H _ , He, He + , and He ++ . Included reactions and their rate coefficients are pre- 
sented in Table 2. 

When the gas density rises, and consequently the chemical reaction timescale becomes shorter 
than the collapse timescale of the cloud, the chemical equilibrium is reached. In addition to this, if 
the cloud becomes optically thick to the continuum, and as a result, the radiation field reduces to 
the black body radiation of the matter temperature, the chemical equilibrium reduce to the Saha 
value. In order to reproduce this feature, inverse processes are included for all reactions. Rates of 
inverse reactions are computed from the principle of detailed balance. 



We switch to equilibrium chemistry at a sufficiently high density n eq where Saha equilibrium 
aeen already reached. In our case, n eq = 10 18 ~ 19 cm~ 3 . 

The rate coefficients of the radiative association 



RA:A + B^C + 7, (CI) 
and its inverse reaction PD, i.e. the photodissociation of C, can be written as 

kKA = K(T)~ X J™ ^^^ ^Dtl - exp(-^)]d,, (C2) 

and 

f°° 4ttJ(u) r . /u/, r <7p n T , /m , 1 n(A)n(B) 11 , 

from the discussion in §A 4.1. 
The partition function of H 

5 

zu = Y,9neM-E n /kT) (C4) 

71=1 
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where g n = 2n 2 and E n = 13.598eV/n 2 . For H2, the sum extends over < v < 5 and < J < 25. 
The rovibrational energies E(v , J) are taken from Borysow et al.(1989). The partition function of 
H2" is taken from Patch k McBride (1968). For other species, we take z e = 2, z H + = 1, z H - = 1, 
z He = 1, z Ue + = 2, and z Uc ++ = 1. 
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Table 1. Continuum Processes 



number 


name 


process 




cross section (cm 2 ) 


al 


H b-f 


Hf n ) + 7 -> H + + e 




7.909 x lCT 18 n(z//z/ n )~ 3 
hv n = 13.598eV/n 2 


a2 


He b-f 


He + 7 -> He + + e 




7.83 x 10- 18 [1.66(^/^ T ) -2 - 05 - 0.66(; 
hv T = 24.586eV 


a3 


H~ b-f 


H" + 7 -> H + e 




10- 18 A 3 (l/A-l/A ) 3 / 2 /(A), A = l 
/(A) given by Eq.(5) of ref. 


a4 


H+ b-f 


H+ + 7 H + H+ 




see Table 2 of ref. 


a5 


H- f-f 


H+e+7^H+e 




kf(T)k B Tn e 

kf(T) given by Eq.(6) of ref. 


a6 


H f-f 


H+ + e + 7^H+ + e 




3.692 x 10 8 g {i (u,T)u- 3 T- 1 / 2 n e , 
we take g&{v,T) = 1 


a7 


H 2 -H 2 CIA 


H 2 («, J)+H 2 + 7^H 2 K 


, ■/') + H 2 


see Fig. 1 of ref. 


a8 


H 2 -He CIA 


H 2 (v, J)+He + 7^H 2 (u / 


, J') + He 


see Fig. 2 of ref. 


si 


H Rayleigh 


h + 7 -» h + y 




5.799 x Kr 29 A" 4 + 1.422 x 10" 30 A- 
+2.784 x KT 32 A~ 8 


s2 


Thomson 


e + 7 — > e + 7' 




6.65 x 1(T 25 



-3.051 



Note. — The wavelength A is in /im. 



References. - (1) Rybicki & Lightman (1979), (2) Osterbrock (1989), (3) John (1988), (4) 
Stancil (1994), (5) Borysow et al. (1997), (6) Kurucz (1970) 
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Table 2. Chemical Reactions 



Number 



Reaction 



Rate Coefficient 



Reference 



1,2 H(n) + e ^ H+ + 2e 



fci >n= i = 4.25 x 10- 11 T 1 /2 exp (- X// /A : T) 

fci,n=2 = 6.69 x 10- 10 T 1 / 2 exp(-Xi//2 2 /^r) (T <4680K) 

fci,n=2: see Janev et al. [2.1.6] (T >4680K) 

fc lin > 3 = 9.56 x 10- 6 T(eV)- L5 

(/3 2 - 33 + 4.38/? 1 / 72 + 1.32/3 n )" 1 exp(-/3„), 
Pn = (xn/n 2 )/kT 
k 2 = kf(z K /z H+ z e )A.UU x 10- 16 T- 1 - 5 exp(x H /A;T) 

K = En=l c n^l,n 



3,4 



H(n) + H ^ H+ + e + H fc ; 



3 — Sn=l c n^3,n 

fc 3) i = 1.2 x 10- 17 T L2 exp(-XH/A;T) 
fe 3 ,n>2 = 7.86 x 10- 15 n 4 /nT a5 (l + 1.27 x 10~ 5 n 2 T) 
(1 + 4.76 x W' 17 n 4 T 2 )- l e^p{-xu/n 2 /kT) 
/ 2 = 0.71, / 3 = 0.81, / 4 = 0.94 
fc 4 = A: 3 *(z h /^h+^)4.1414 x 10- 16 T- 1 - 5 exp( X H/A ; T) 

^3 = En=l c n^3,n 



1 

3 



5,6 



H- + H+^H(n)+H 



^5 — En=l ^5,n 

fc 5) i = 6.92 x 10' 14 T - 5 
% 2 = 8.0 x 10- 13 T - 83 
fc 5 ' 3 = 6.18 x 10- 7 T-°- 27 

h = [J2l =1 (Cn/c* n )h, n ] 

(z H -z H+ /z H 2 )exp[-(xH - Xn-)/kT] 



7,8 



H 2 + H+ ^ H + 



k 7 = 1.5 x 10- 10 exp(-1.4 x 10 4 /T) (T > 10 4 K) 
fc 7 = 3 x 10~ 10 exp[-( X H 2 - Xu 2 +)/kT] (T < 10 4 K) 
k 8 = /c 7 (z H+ ZH 2 AH% 2 +)exp[(xH 2 - Xu 2 +)/kT] 



9,10 



H + H" 



H 2 + e 



fc 9 = 1.5 x 10~ 9 (T < 3 x 10 2 K) 

fc 9 = 4.0 x l(r 9 T- a17 (T > 3 x 10 2 K) 



k 



"-in tr"'i n 



10 — n '10,H n '10,L 

fci 0)L = 2.7 x lO-ST-^exphCxHa 



x H -)An 
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Table 2 — Continued 



Number Reaction Rate Coefficient Reference 

fcio,H = k 9 (znz u -/zu 2 z e )2.775 x 10 4 exp[-(xH 2 ~XH-)/kT] 
a = (1 + n/n cr ) _1 

log 10 (n cr ) = 4.0 - 0.4161og 10 (T/10 4 ) 5 
-0.327(log 10 (T/10 4 )) 2 

11,12 H+ + e^H(n)+H k n = Y? n =i k n , n = 2 x lO^T' 1 / 2 4 

ku,i : fcii )2 : fcii,3 = &11.4 : = : 0.10 : 0.45 : 0.22 : 0.12 2 
but, k 11:n = in case of xh - X H + ~ Xu/n 2 < kT 

kl2 = En=l( c n/ c n) fc H,n] 

( %2+ z e /z H 2 )3.6034 x 10- 5 exp[-( XH -Xn 2 +)/kT] 

13,14 2H 2 ^ 2H + H 2 k 13 = fc^^is.L 5 

jfe 13)L = 1.18 x 10- 10 exp(-6.95 x 10 4 /T) 
fc 13)H = 1.30 x 10- 9 exp(-5.33 x 10 4 /T) 
a = (1 + n/n CT )~ l 

log 10 (n cr ) = 4.845 - 1.31og 10 (T/10 4 ) 
+ 1.62(log 10 (T/10 4 )) 2 
ku = £:i3,h(W^h 2 )1-493 x 10- 20 T- L5 exp(xH 2 AT) 

15,16 H 2 + H ^ 3H k 15 = k\^kf 5jL 

fci 5 ,L = 1-12 x 10- 10 exp(-7.035 x 10 4 /T) 5 
jfe 15iH = 6.5 x lO -7 ! 1-1 / 2 6 

exp(-XH 2 /^T)[l - exp(-6000/T)] 
a = (1 + n/n cr ) _1 

log 10 (n cr ) = 4.0 - 0.4161og 10 (T/10 4 ) 5 
-0.327(log 10 (T/10 4 )) 2 
fcie = A: 15 ,h(^h 2 Ah 2 )1.493 x 10- 20 T- 1 - 5 exp( X H 2 /A;T) 

17,18 H 2 + e ^ 2H + e k 17 = 1.3 x 10- 18 T 2 exp(-XH 2 A^) 1 

fcis = MWV) 1.493 x lO-^T-^exp^AT) 

19,20 H- + e ^ H + 2e fc 19 = 4 x 10- 12 Texp(-x H - 5 

k 20 = k 19 (z n -/z R z e )4.UU x 10- 16 r- 1 - 5 exp( XH -/A;T) 
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Fig. 1. — Two FUV spectra studied in this work. The solid and dotted lines indicate those of Types 
a and b, respectively. The normalization of intensity J 2 i = 1 for both types. Also shown are the 
threshold frequency of photodissociation (0.755 eV) and a representative value of Lyman and 
Werner (LW) bands (12.4 eV). 

Fig. 2. — The density-temperature relation for the collapse of primordial clouds irradiated with FUV 
radiation. The spectra are type a: Juv(^) = J21 x 10~ 21 (v / Vth) -1 ■ The evolutionary trajectories 
are drawn for J 21 = 0, 10~ 2 , 1, 10 2 , 10 4 , and 10 5 . The trajectories for J 21 > 10 5 are identical to 
that for J21 = 10 5 . The dashed lines indicate the constant Jeans mass. The dark matter gravity is 
neglected in calculating the Jeans masses. 

Fig. 3. — Same as Fig. 2, but for spectra type b: J\jv( v ) = J21 x 
10~ 21 [B(i/; W 4 K)/B(u t i l ; 10 4 K)] [y < u^). The evolutionary trajectories are drawn for J 2 i = 
0, 10 -2 , 1, 10 2 , and 10 3 . The trajectories for J 21 > 10 3 are identical to that for J 21 = 10 3 . 

Fig. 4. — Concentration of (a) hydrogen molecules y(H 2 ) and (b) electrons y(e) for the clouds 
irradiated with the power-law type radiation (type a spectrum in the text) with J21 = 10 5 (solid 
lines), 10 4 (dotted lines), 10 2 (short- dashed lines), 1 (long-dashed lines), 10~ 2 (dash-dotted lines), 
and (long- dash- dotted lines). 

Fig. 5. — Cooling/heating rates per unit baryonic mass as a function of the central number density 
for the clouds irradiated with the power-law type radiation of (a) J 2 i = 0, (b) 10 4 and (c) 10 5 . 
These include the contributions by the compression, chemical reactions, atomic hydrogen lines, H 
two-photon emission, continuum, and H 2 lines. The dominant continuum processes are (a) H~ 
free-bound emission (n < 10 11 ' 5 cm~ 3 ), H 2 collision-induced emission (10 11-5 cm -3 < n); (b) H~ 
free-bound emission (n < 10 11 ' 5 cm~ 3 ), H 2 collision-induced emission (10 1L5 cm -3 < n); and (c) H~ 
free-free emission (10 2 < n < 10 4 ' 5 cm -3 ), H~ free-bound emission (10 45 cm -3 < n), respectively. 
The chemical reactions contributing dominantly to the cooling/heating rates are: (a) and (b): 
H 2 formation (n < 10 13 cm~ 3 ), H 2 dissociation (10 13 < n < 10 15 cm -3 ), H 2 formation (10 15 < 
n < 10 16 cm -3 ), H 2 dissociation (10 16 < n < 10 21 cm -3 ), H ionization (10 21 cm~ 3 < n), and 
(c) H 2 formation (n < 10 155 cm" 3 ), H 2 dissociation (10 15 - 5 < n < 10 16 - 5 cm~ 3 ), H ionization 
(10 16 ' 5 cm -3 < n), respectively (cf. Fig. 4). 

Fig. 6. — Estimated H 2 fraction vs. the amount of H 2 needed for efficient H 2 cooling during the 
initial adiabatic temperature rise. We assume that the temperature rises adiabatically from the 
initial state of calculation, that is, T = To(n/no) 2 ^ 3 , where To = 39 K, and no = 8.9 x 10~ 2 cm~ 3 
are the initial temperature and number density, respectively. The ionization degree is taken to be 
10~ 4 . The H 2 fractions are estimated by equation (A2). The solid and dotted lines illustrate the 
H 2 fractions for Types a and b FUV radiation whose values of J 2 i are indicated in the figure. The 
thick solid curve shows the necessary H 2 fraction for the cloud to cool in a free-fall time (eq. [Al]) 

Fig. 7. — Estimated H 2 fraction vs. the amount of H 2 needed for efficient H 2 cooling for clouds 
that collapse isothermally by atomic cooling. The temperature and ionization degree are taken to 
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Table 2 — Continued 



Number Reaction Rate Coefficient Reference 

21,22 + e ^ H + H+ + e k 2 i = 2 x l(T 7 exp(-x H + /kT) 1 

k 2 2 = k 2 i(z H +/z a z H+ ) 1.493 x 10- 20 r- L5 exp(x H +/feT) 

23,24 He + e ^ He+ + 2e fc 23 = 2.38 x lO'^T 1 / 2 ^ + (T/10 5 )]- 1 exp(-XHc/A:T) 7 

fe 24 = MW%e+*e) 4 - 1414 x lO-^T'^eMXHe/kT) 

25,26 He+ + e ^ He++ + 2e fc 25 = 5.68 x Hr^T 1 / 2 ^ + (T/10 5 )]- 1 exp(-x H c+/^) 7 

fee = fe5(% c +/% c ++^)4.1414 x 10- 16 r- 1 - 5 exp(x Hc +/feT) 

27,28 H + 7 ^ H+ + e see equations (C2), (C3) in the text 

and process (al) in Table 1 

29,30 He + 7 ^ He + + e see equations (C2), (C3) in the text 

and process (a2) in Table 1 

31,32 H~ + 7 ^ H + e see equations (C2), (C3) in the text 

and process (a3) in Table 1 

33,34 H+ + 7 ^ H + H+ see equations (C2), (C3) in the text 

and process (a4) in Table 1 

35 H 2 + 7 -► W 2 ^ 2H k 35 = 1.4 x 10 9 J(hv = 12.4eV)/ sh 8 

/ sh = min[l,(iVH 2 /10 14 cm- 2 )- 3 / 4 ] 



Note. — The temperature T is in K, except otherwise noted. 

Note. - The binding energies are xn/k = 1.578 x 10 5 K, xh 2 A = 5.197 x 10 4 K, Xh-/& 
.761 x 10 3 K, Xll +/k = 3.067 x 10 4 K, X Hc/k = 2.853 x 10 5 K, and xhc+A = 6.312 x 10 5 K. 



References. — (1) Lenzuni et al. (1991), (2) Janev et al. (1987), (3) Drawin (1969), (4) Galli & 
Palla (1998), (5) Shapiro & Kang (1987), (6) Palla et al. (1983), (7) Black (1981), (8) Draine & 
Bertoldi (1996) 
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be 8000K and 10~ 4 , respectively. The H2 fractions are estimated by eq.(A2). The solid and dotted 
lines illustrate the H2 fractions for types a and b FUV radiation whose values of J21 are indicated in 
the figure. The thick solid curve shows the necessary H2 fraction for the cloud to cool in a free-fall 
time (eq.[Al]). 
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